Analysis of the Superkamiokande atmospheric neutrino data in the 
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Superkamiokande atmospheric neutrino data for 990 days are analyzed in the framework of four neutrinos 
without imposing constraints of Big Bang Nucleosynthesis. It is shown that the wide range of the oscillation 
parameters is allowed at 90% confidence level (0.1 <|(7 s i| 2 + \U S 2\ 2 < 1). 



1. Introduction 

It has been known that three different kinds of 
experiments suggest neutrino oscillations: the so- 
lar neutrino deficit fij the atmospheric neutrino 
anomaly ||[| and the LSND data (§. If one as- 
sumes that all these three are caused by neutrino 
oscillations then one needs at least four species 
of neutrinos. Furtherover, it has been shown 
(5Hl that the 4 x 4 MNS matrix splits approx- 
imately into two 2x2 block diagonal matrices if 
one imposes Q] the constraint of the reactor data 
JsJ and if one demands that the number N v of 
effective neutrinos in Big Bang Nucleosynthesis 
(BBN) be less than four. In this case the solar 
neutrino deficit is explained by v e <-> v s oscilla- 
tions with the Small Mixing Angle (SMA) MSW 
solution and the atmospheric neutrino anomaly is 
accounted for by «-> v T . On the other hand, 
some people have given conservative estimate for 
N v H and if their estimate is correct then the 
constraints on the mixing angles of sterile neutri- 
nos are not as strongs as in the case of N v < 4 
and the conditions one has to take into account 
are the data of the reactor, the solar neutrinos 
and the atmospheric neutrinos. Recently Giunti, 
Gonzalez-Garcia and Peha-Garay JnJ have an- 
alyzed the solar neutrino data in the four neu- 
trino scheme without BBN constraints. They 
have shown that the scheme is reduced to the two 
neutrino framework in which only one free param- 



eter c s = \U s i\ 2 + \U S 2\ 2 appears in the analysis. 
Their conclusion is that the SMA MSW solution 
exists for the entire region of < c s < 1, while the 
Large Mixing Angle (LMA) and Vacuum Oscilla- 
tion (VO) solutions survive only for < c s <0.2 
and < c s <0.4, respectively. In this talk I will 
discuss the analysis of the Superkamiokande at- 
mospheric neutrino data for 990 days [|| (con- 
tained and upward going through /x events) in 
the same scheme as in ]To[ |, i.e., in the four neu- 
trino scheme with all the constraints of acceler- 
ators and reactors but without BBN constraints. 
Details and more references are given in |l]J| . 



2. Four neutrino scheme 

Here I adopt the notation in 
MNS matrix: 

U e2 U e3 
U T2 U T3 

U S 2 U s3 



for the 4 x 4 
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I can assume without loss of generality that 
m\ < m\ < m§ < m|. Three mass scales 
Am% ~ O(10- 5 eV 2 ) or 0(lO~ lo eV 2 ), Am 2 tm 
- 0(lO- 2 eV 2 ), Am 2 SND ~ C(feV 2 ) are neces- 
sary to explain the data of the solar neutrinos, 
the atmospheric neutrinos and LSND, so I as- 
sume that three independent mass squared dif- 
ferences are Am^, Am 2 tm , AmL SND . It has 
been known (|,[l^j that schemes with three de- 



generate masses and one distinct massive state 
do not work to account for all the three neu- 
trino anomalies, but schemes with two degen- 
erate massive states (m 2 ~ m 2 <§; m§ ~ m 2 , 
where (a) (Amfi, Aml 3 ) = (Atoq, Am 2 tm ) or 
(b) (Am 2 21 ,Amj 3 ) = (Am^, Am|)) do. As far 
as the analyses of atmospheric neutrinos and solar 
neutrinos are concerned, the two cases (a) and (b) 
can be treated in the same manner, so I assume 
Am^ = Attjq, Am 2 3 — Am^ in the following 
(See Figg). 
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Figure 1. Four neutrino scheme 



For the range of the Am 2 suggested by 
the LSND data, which is given by 0.2 
eV 2 <Am 2 SND < 2 eV 2 when combined with the 
data of Bugey JsJ and E776 Jl3| , the constraint 
by the Bugey data is very stringent and 



\U e , 



|iU 2 <io- 2 



(2) 



has to be satisfied ]7|j5|,^2||. Therefore I put 
U e3 = [/ e 4 = for simplicity in the following 
discussions. Also in the analysis of atmospheric 
neutrinos, \Am 2 7) L/AE\ <C 1 is satisfied for typi- 
cal values of the neutrino path length L and the 
neutrino energy E, so I assume Am|x = for 
simplicity throughout this talk. 

Having assumed C/ e 3 — U e ^ — and Am^ = 0, 
I have only mixings among , v T , v a in the analy- 
sis of atmospheric neutrinos, and the Schrodinger 



equation I have to consider is 

d ( ""^ 
X V v s {x) 



M 



t/diag(-A J B 3 2,0,A£43) 
+diag(0,0,A(^))] 



u T (x) 
v s {x) 

u- 1 



(3) 



where AE^ = Amf j /2E, A(x) = G F N n {x)/^/2 
stands for the effect due to the neutral current in- 
teractions between v T and matter in the Earth 
and 



U = 




D e l 



(4) 



with D ee diag (e iSl / 2 , 1, e"^ 1 / 2 ) (Aj are the 3 x 3 
Gell-Mann matrices) is the reduced 3x3 MNS 
matrix. This MNS matrix U is obtained by sub- 
stitution 6> 12 — > 6>23 - 7T/2, 9 13 — > 6*24, 6*12 ~> 

7r/2 — #34, 5 — ■+ Si in the standard parametriza- 
tion in [JlS. Since v e does not oscillate with any 
other neutrinos, the only oscillation probability 
which is required in the analysis of atmospheric 
neutrinos is P(v^ — > v^). 

In the following analysis I will consider the situ- 
ation where non-negligible contribution from the 
largest mass squared difference Am\ 2 appears in 
the oscillation probability P{t>^ — * v^). In short 
baseline experiments where |Ai?2i| and IA-E43I 
can be neglected, the disappearing probability 
P(^ M ^ is given by 

- = l-4|t^ 2 | 2 (l-|^2| 2 ) 



ATO32L 

AE 



(5) 



so that the mixing angle sin" 2#sbl is given by 

sin 2 20 S BL = 4c 2 4S 2 3 (1 - c 2 4S 2 3 ) . (6) 

It turns out in the final results that sin 2 2#sbl 
can be as large as 0.5 in the allowed region of 
the atmospheric neutrino data. To avoid contra- 
diction with the negative result of the CDHSW 
disappearing experiment on — > |l5| ], I will 
take Am 2 2 =0.3eV 2 as a reference value. 
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Figure 2. Allowed region at 90%CL for 5 X = tt/2. 

The shadowed area is the allowed region projected on the (#34, #23) plane for various values of A 771,43 
(10~ 3 5 eV 2 < Am| 3 < 10~ 2 cV 2 ) for each value of # 24 = 30°, • • • , 55°, respectively, and the thin solid lines 
are boundary of the allowed region for various values of Ara 2 3 . The solid, dashed, coarse dotted and fine 
dotted lines stand for the contours of c s = \U s i\ 2 + \U S 2\ 2 HC23C34 + S23S24S 3 4 : e lSl | 2 = 0.2,0.4,0.6,0.8, 
respectively. Solutions with c s <0.2 exist for 40° <^24 <55° and they can have Large Mixing Angle 
solutions of the solar neutrino problem. 



3. Analysis of the atmospheric neutrino 
data 

I calculate the disappearance probability 
P{ v fj, — > ^fj.) by solving (^) numerically, and eval- 
uate the number of events. Then I define x 2 a s 



2 _ 2 

X — Xsub-GeV 



Xmulti — GeV ' Xthrough (7) 
- ^2 Qr ,j „2 



U = -6,---,6) 
5i = 0°,45 o ,90° 



wh ere X sub _ GeV , X mu iti-GcV and Xthrough are 
% 2 for sub-GeV, multi-GeV, and upward going 
through /i events, respectively. I have evaluated 
X 2 for 6 2i = (25 + 5j)° (i = 0, • • • , 7), 6- 3 4 = 15j° 
3 = 10j° (j = 0,---,4), 
m\ 3 = 10- 4 +^ 10 eV 2 (j = 
5, • • • , 20) and it is found that x 2 nas the min- 
imum value Xmin = 43- 1 (Xsub-GcV = 19 A 

Xmulti-GeV = 13 ' 2 > Xthrough = U - 6 ) for Am 43 = 
l -2.9 eV 2 = 13 x 10 -3 eV 2 j (e 2 4,9 34 ,9 23 ) = 

(35°, 15°, 20°), Si = with 45 degrees of free- 
dom. For pure «-> v T (634 — 623 — 0), 
the best fit is obtained Xmin(^M ^ v t) = 48.3, 

(Xsub-GcV = 19.8, Xmulti-GcV = 17.0, Xthrough = 

10.6) for Aml 3 = 2.0 x 10~ 3 eV 2 , (0 24 , 6» 3 4, M = 
(40°,0°,0°). The allowed regions at 90%CL are 
obtained by x 2 < Xmin + A X 2 : where Ax 2 =9.2 for 
five degrees of freedom. In Fig.^ the allowed re- 
gion at 90% confidence level is depicted as a shad- 
owed area in the ($34, #23) plane for various val- 
ues of Aml 3 (10~ 35 eV 2 < Aml 3 < 10~ 2 eV 2 ) for 
each value of #24 — 30°, • • • , 55° and for Si = tt/2, 
together with lines c s =constant. A few remarks 
are in order. (1) Pure <-> v s oscillation, which 
is given by #34 = ±90°, $23 = 0, is excluded 
at 99.7%CL for any value of Am 2 ,, #24, Si and 



this is consistent with the claim ]3j by the Su- 
perkamiokande group. (2) For generic value of 
($34, $23), the oscillation is hybrid not only with 



v T and Vu 



but also with Ar 



'43 



and 



Am 32 . To illustrate this, let me give the expres- 
sions of oscillation probability in vacuum in the 
case of Si = tt/2: 



„2 „2 ,2 „2 
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where I have averaged over rapid oscillations due 
to Am 32 - As is seen in (||), roughly speak- 
ing, #34 represents the ratio of <-> v T and 
Vft «-> v Sl whereas 6*23 indicates the contribution 
of sin 2 (Am 32 L / AE) in oscillations. Zenith angle 
dependence of the /x-like multi-GeV events and 
the upward going through \i events are shown 
in FigJ^ for a few sets of the oscillation parame- 
ters. The disappearance probability behaves like 
1 - P(z^ -»■ Vp) = a + [3 sin 2 (Am 2 i3 L / AE) (a, 
f3 are constant) and a > is satisfied when- 
ever #23 7^ 0. Because of this constant a, which 
never appears in the analysis of the two flavor 
framework, the fit for 6*23 ^ tends to be bet- 
ter than in the case of #23 = 0. The reason 
why the best fit point is slightly away from pure 
<-» v T case and the reason why an exotic solu- 
tion like (6> 24 , 6*23, 3i ) =(45°, 30°, 90°), S x = 90°, 
Am| 3 = 1.3x 10 _3 eV 2 is allowed is because a bet- 
ter fit to the multi-GeV contained events com- 
pensates a worse fit to the upward going through 
\x events, and in total the case of hybrid oscilla- 
tions fits better to the data (Notice that the fit 
of <-> v s scenario to the contained events is 
known to be good |l6|,[l7j and in the present case 
the fit becomes even better due to the presence 
of a). 



4. Conclusions 

I have shown in the framework of four neu- 
trino oscillations without assuming the BBN con- 
straints that the Superkamiokande atmospheric 
neutrino data are explained by wide range of the 
oscillation parameters which implies hybrid os- 
cillations with «-> v r and «-> v s as well 
as with Am 2 tm and Am 2 SND . The case of pure 
Vfi v s is excluded at 3.0<rCL in good agreement 
with the Superkamiokande analysis. It is found 
by combining the analysis on the solar neutrino 
data by Giunti, Gonzalez-Garcia and Peha-Garay 
that the LMA and VO solutions as well as SMA 
solution of the solar neutrino problem are allowed. 
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Figure 3. Zenith angle dependence 
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